Sling effect in collisions of water droplets in turbulent clouds 
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We describe and evaluate the contribution of sling effect into the collision rate of the same-size 
water droplets in turbulent clouds. We show that already for Stokes numbers exceeding 0.2 the sling 
effect gives a contribution comparable to Saffman- Turner contribution, which may explain why the 
latter consistently underestimates collision rate (even with the account of preferential concentration). 
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PACS numbers: 

The collision rate N is the product of the target area 
Tr(a+a') 2 , the relative velocity of droplets before the con- 
tact Av and the probability P to find two droplets touch- 
ing: 



N(a, a') = ir(a+a') 2 AvP{a + a') 



(1) 



where a, a' are the radii of the droplets. Here we focus 
on the contributions of the air flow to the collision rate. 
To that end we consider equal-size droplets which fall 
with the same velocity in still airj2^| We also consider 
droplet sizes exceeding few microns and neglect Brownian 
motion of droplets. For such droplets, the air flow is the 
sole source of relative velocity Av and it also influences 
P(2a) due to droplet inertia. 

Let us first briefly discuss the latter effect called pref- 
erential concentration. It has been identified long ago 
(see Maxey 1987; Squires and Eaton 1991; Sundaram 
and Collins 1997; Reade and Collins 2000; Kostinski 
Shaw 2001; Jaczewski and Malinowski 2005; McFarquhar 
2004; Franklin et al 2005; Grits et al 2006 and the ref- 
erences therein). Still, the proper quantification of this 
effect and of its role in the collision rate enhancement 
in warm clouds remains to be done. It has been in- 
ferred from the data (Sundaram and Collins, 1997) that 
P(l) has a power-law dependence and argued theoreti- 
cally (Balkovsky et al, 2001) that the dependence must 
be of the form P(l) ~ {r]/l) a where r\ is of the order 
of the viscous scale of turbulence and the dimensionless 
quantity a depends on the dimensionless numbers that 
characterize air turbulence, particle inertia and gravity 
— Reynolds, Stokes and Froude numbers respectively: 

Re = (£/??) 4/3 , St = Tv/rj 2 = At , F = Xrj/gr . (2) 

Here L is the size of the largest turbulent eddies (typ- 
ically, the size of the cloud), v is the air viscosity, g is 
the acceleration of gravity and r= (2/9)(po/ p)(a 2 /v) is 
called the reaction (Stokes) time with p, po being the air 
and water density respectively. 

The behavior of a is well-understood for small St where 
a ~ b(Re, F)St 2 as predicted by Balkovsky et al (2001), 
Falkovich et al (2002), and confirmed by Falkovich and 
Pumir (2004) and Chu and Koch (2005) for moderate 
Re. Here we establish a for arbitrary St and substan- 
tially higher Re than before. Note in passing that a gen- 
eralization for different-size droplets value, P{a,a'), has 



been suggested in Falkovich et al (2002) and Bee et al 
(2005). 

The main subject of the paper is the proper evalua- 
tion of the contribution of the relative velocity into the 
collision rate, particularly, the (generally nonlocal) re- 
lation between the air flow and the droplet velocities. 
Saffman and Turner (1956) assumed that the relative ve- 
locity of the droplets is determined locally by the air 
velocity which is spatially smooth at such scales (since 
a ^ rj) so that Av ~ 2Aa where A is the rms air-velocity 
gradient according to (|2Jl. However, the droplet velocity 
is determined not only by a local air velocity but also 
by the previous history because of inertia. As was first 
noticed by Falkovich et al (2002), this leads to an ex- 
tra contribution to the collision rate not captured by the 
Saffman- Turner formula. Figure 1 illustrates this (so- 
called sling) effect: The right droplet passed through 
an intense vortex and had been thrown away as if by 
a sling. As a result, the relative velocity of the droplets 
at the point of collision may be determined not by the 
air-flow gradient at this point but rather by a distant 
vortex. Let us give here some numbers adopted from 
the calculations to be described below. Consider turbu- 
lence with A = 80s -1 . For droplets with a — 10pm, 
St = 0.08 and sling effects are negligible; for a = 15 pm, 
St = 0.2 and the frequency of sling effects is ~ A/1000 
while their contribution into the collision rate is about 
20%; for a = 20 pm, St = 0.35 and the frequency of 
sling effects is ~ A/50 while their contribution into the 
collision rate is about 35% (the numbers correspond to 
Figs. 3,5,6 below). An independent confirmation that 
Saffman- Turner formula (which disregards sling events) 
consistently underestimates the collision rate even at rel- 
atively small Stokes numbers has been obtained recently 
by direct numerical simulations (Franklin et al, 2005). 

The most straightforward way to model the collision 
rate is with the discrete finite-size particles embedded 
into a properly modelled flow (see, for instance Franklin 
et al, 2005 and Wang et al, 2005). Since turbulent air 
flow and droplet distribution are highly intermittent, the 
requirements on resolution and statistics are such that 
both Reynolds numbers and number of droplets are mod- 
erate at best in such a modelling. Moreover, compu- 
tations with discrete particles do not allow one to dis- 
tinguish between local contributions to collisions (that 
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FIG. 1: Sketch of the sling effect. Dotted lines show the tra- 
jectories of the droplets while broken lines show the stream- 
lines of the vortex in the air flow. 

must be described by the Saffman- Turner formula with a 
proper preferential concentration correction) and a sling- 
effect contribution. That impedes the progress towards 
a proper parametrization of the collision rate. Here we 
propose another, complementary way of modelling based 
on the continuous description of the flow of droplets. The 
equation for the droplet velocity v, 

dv/dt = (u - v)/r + g , (3) 

can be considered as defining everywhere in space the 
field v(r) from the known air-flow field u(r). In the con- 
tinuous description we are able to see clearly the sling 
contribution. Indeed, sling events appear in the equa- 
tions of motion for particles as crossing of trajectories or 
finite-time singularities of the continuous equations. For 
the matrix of droplet-flow gradients, er^ = dvi/drj, one 
gets in the co-moving (Lagrangian) reference frame: 

Vij + 0ik&kj + "Xj/t = Sij/r , (4) 

where Sij = dui/drj. Smoothness of the air flow means 
finiteness of Sij(t). On the contrary, nonlinearity (which 
corresponds to inertia) in the equation for a leads to 
the possibility of explosions when some component of 
turns into -oo in a finite time by the law a cx (t — to) -1 
(Falkovich et al, 2002). 

Because of spatial and temporal non-locality, ana- 
lytical description of the sling-effect contribution into 
the collision rate is difficult. Using simple models, it 
has been predicted that probabilities of such events 
have a very sharp dependence on the Stokes number: 
exp[-c(Re,F)/St] (Falkovich et al 2002; Wilkinson and 
Mehlig 2003,2005), where neither the value of the factor 
c(Re, F) nor its dependence on the Reynolds and Froude 
numbers is known. Particularly important is understand- 
ing the dependence on Reynolds number which varies by 
many orders of magnitude in atmospheric flows. 

In this paper, we perform direct numerical simulations 
of the air-flow turbulence at such high Re that were never 
reached before in collision rate calculations. We pay the 



price by not performing kinetic droplet simulations (i.e. 
following separate droplets) but by using continuous field 
description based on (|3I4|I and the equation for particle 
concentration in the Lagrangian reference frame, 

dn/dt = —nan(t) . (5) 

To this end, we generate a statistically stationary turbu- 
lent flow in a cube with periodic boundary conditions. 
The equations are solved with a pseudo spectral code, 
see Pumir (1994) for details. All the appropriate length 
scales of the flow are adequately resolved. We work in the 
range 21 < R\ < 105. In this turbulent flow, we follow 
the motion of incrtial particles by solving the equation 
for the position, dx/dt = v, along with ©. The equa- 
tion 0| for the tensor of droplet velocity derivatives is 
integrated along the way. The integration of Eq. 13141 
requires the interpolation of the fluid velocity, u, or its 
derivative, s, from the numerical mesh to the particle 
position. This is done by using spline interpolation tech- 
niques, see Girimagi and Pope (1990). The equations are 
solved by using algorithms that are second order accurate 
in time, or higher. 

For each run, we integrate the equations of motion 
for several Stokes numbers (ranging from St = 0.05 to 
St « 5). Gravity is also taken into account, by taking 
a finite value of the Froude number. In this work, as 
in Falkovich and Pumir (2004), the only property of the 
droplets which varies is the radius, a. Accordingly, the 
Froude and the Stokes number vary in such a way that 
the product of the Stokes number by the Froude number 
is constant : 

e = St x Fr (6) 

Numerically, it was found that eo ^ 5 corresponds to a 
vanishingly weak gravity. The effect of gravity becomes 
appreciable for values of eo smaller than 1. 

The occurrence of a sling effect results in a singularity 
of a in a finite time, a cx (t — to) , which in turn, leads 
to a divergence of the particle density (n cx (t — to) -1 )- 
Physically, neither the droplet velocity gradient nor the 
droplet density can grow unrestricted since droplets have 
a finite size, a, and cannot come arbitrarily close to one 
another. Therefore, the droplet velocity gradients do not 
grow by more than by a factor ~ (rj/a). Once the gradi- 
ent has reached this predetermined threshold value, the 
equation is regularized, by flipping the sign of a. This 
simply corresponds to a fast droplet passing a slow one so 
that their velocity difference changes sign (as well as the 
velocity gradient). Similarly, the density increases until 
the time of the flip and then decreases. This algorithm 
leads to a numerically well-posed problem. 

To infer the properties of the coarse grained distribu- 
tion at a given scale r, we use the method developed by 
Balkovsky et al (2001), Falkovich et al (2002) and imple- 
mented numerically by Falkovich and Pumir (2004) in the 
restricted case where the Stokes number was small, and 
where no sling contribution was expected. That method 
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requires one to follow (in addition to x, v, a and n) the 
deformation of a volume, seeded with bubbles, and car- 
ried by the v-flow. Specifically, one needs to determine 
the contraction rate of the volume along particle trajec- 
tories. This is effectively done by monitoring the growth 
of the inverse of the deformation tensor, W, which de- 
scribes how a line element is transported by the flow : 
5l{t) = W(t) ■ (51(0), where W' 1 satisfies : 

( M^L = -( W - 1 -a + a T -W- 1 ) (7) 

To estimate the contribution of a trajectory to the coarse 
grained particle density at scale r, the integration is car- 
ried out until W^ 1 reaches (r)/r). At this point, the value 
of n is recorded. The Eulerian value of the k th moment 
is obtained by averaging the values of n k ~ 1 over all the 
trajectories computed. 

In such an approach, the collision rate along a trajec- 
tory consists of two contributions. The first one is the 
'continuous contribution', determined by a local velocity 
gradient (like in the Saffman- Turner approach). Given 
the value of n and of a at a given time, the instanta- 
neous flux of incoming particles towards a given particle 
is as follows: 

$cont(t) = -(2afn(t) I (e-a-e)dfi (8) 

Jecre<0 

The contribution to the collision rate, K cont , along a tra- 
jectory is simply obtained by integrating $ con t over time. 
Because the growth of n is accounted for in the result- 
ing contribution for the collision term, the influence of 
the sling effect on concentration (because of caustics left 
after sling events) is taken into account in this formula. 
However, as has been pointed out above, the sling effect 
also results in an additional contribution to the velocity 
difference, which is not proportional to the local velocity 
gradient. As suggested by Fig.l, a sling event involves 
a number of particles incoming in a small region with 
significantly different velocities, a situation leading to an 
outbursts of collisions, which we estimate as follows. The 
source term in Eq. 0] necessary to start a blow-up pro- 
cess, s/t, should be large enough to drive a to start the 
blow up : |s| > 1/r. Based on the Kolmogorov picture of 
fully developed turbulence, the extent of the region where 
gradients reach the value s of the order or larger than 1/r 
is of the order of I ~ {vt) 1 !" 1 = -qSt 1 / 2 . The time over 
which the collision takes place is estimated to be of the 
order of the droplet relaxation time, r. Last, the range 
of particle velocities involved during the collision can be 
estimated as \8w\ ~ 1/t. Based on these estimates, the 
number of collisions that occur in the wake of a sling 
event that have happened at time t s can be estimated as 
follows: 

N sUng (t) = 47r(2a) 2 x n(t s + r/2) x |<5v| x r (9) 

The relaxation time of droplet velocities is r and this is 
a typical duration of the time when droplets have sub- 
stantially different velocities and sling-effect contribution 



appears. That is why in estimating N s u ng the density of 
particles is taken at a time r/2 after the sling event time 
t s , which provides a reasonable estimate for the particle 
density during the entire process. As it is clear from the 
derivation, the value of N s u ng obtained in such a way 
is up to a numerical factor of order unity, whose precise 
value could be estimated by kinetic numerical simula- 
tions, which is beyond the scope of this work. In the rest 
of this paper, we estimate separately the contributions 
from the regular term and the sling contribution, given 
by Eq. |S] and Eq. respectively. 

The method sketched above to estimate the coarse 
grained properties at a scale r is then used to evaluate 
the collision rate as a function of scale. More precisely, 
consider droplets of radius a. To compute the collision 
rate, the regular and sling contributions to the collision 
rate are computed along trajectories, until the compres- 
sion, given by see Eq.[7| reach the value i)]a. The 
various contributions coming from different trajectories 
are then accumulated, and the mean value of the collision 
rate extracted. 

The way the flow leads to the compression of an en- 
semble of particles, described by Eq. [7| plays a crucial 
importance in the physical processes controlling parti- 
cle collision rates. The numerical results indicate that 
ZndPF -1 )) grows linearly with time. The growth rate 
gives a direct access to the smallest Lyapunov exponent 
A3 that corresponds to the strongest contraction in the 
flow. At a given Reynolds number, the value of —A3 x tk 
has a non trivial dependence on the Stokes number. It 
starts at a value —A3 x tk ~ 0.16 at St = (Falkovich 
et al, 2002), then increases up to St ~ 0.5, leading to a 
twofold increase of A3 compared to its value when St = 0, 
before decreasing at higher values of St, see Fig. 2. Qual- 
itatively similar dependence has been found for a model 
short-correlated flow (Wilkinson and Mehlig 2003). The 
values of A3 are not affected much by the effect of gravity, 
at the moderate values of the Froude numbers considered 
here, as shown in the figure. No significant dependence 
of the product tk x A3 as a function of the Reynolds 
number was found over the range covered here. 

Sling events, in our approach, are manifested by diver- 
gences of the particle velocity derivative tensor, a. The 
blow-up frequency, fbu, defined as the total number of 
sling events, divided by the integration time, plays here 
a crucial role. Fig. 3 shows fbu multiplied by the Kol- 
mogorov time tk- The data shown correspond to two 
values of eq (eq = 5 and eo = 0.4). The blow-up fre- 
quency, made dimensionless with the Kolmogorov time, 
tk, shows a fairly weak dependence of the Reynolds num- 
ber. No blow-up is observed at very low value of St (for 
St < 0.15). The value of fbu then raises to a maxi- 
mum value for St ~ 1.5, then decreases slowly. Upon 
increasing gravity (decreasing eo), the blow-up frequency 
generally goes down, with a similar Stokes number de- 
pendence. The recent works on simple versions of the 
problem (1-dimensional and short-correlated flows sug- 
gest the dependence of the blow-up frequency as a func- 
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FIG. 2: The contraction rate along trajectories of particles 
at R\ = 105 and at a very low gravity (eo = 5, upper curve) 
and at a moderate gravity (eo = 0.4, lower curve). 

tion of St of the form: fi, u » exp(— A/ St) (Wilkinson 
and Mehlig 2005, Derevyanko et al 2006). Here we find 
empirically that the curve could be fit pretty well by the 
dependence of the form: 

fbu y-T K = Sr 2 x exp(~A/St) x (B + CSt c ) (10) 

The coefficient A is found to decrease slightly as the 
Reynolds number increases (A — 2.1 for R\ = 45, 
A = 1.85 for R x = 83 and A = 1.70 for R x = 105), con- 
sistent with the fact that as turbulence becomes more 
intense, higher gradients appear in the flow, which are 
able to induce blow-up of a at increasingly low values 
of St. Yet, at higher values of the Stokes number, the 
blow-up frequency seems to decrease as the value of the 
Reynolds number increases, a somewhat surprising effect. 

The dependence of the coarse-grained particle density, 
(n 2 ) r as a function of rj/r, is very similar to the one 
obtained by Falkovich and Pumir, 2004. Namely, (n?) r 
has essentially a power-law dependence as a function of 
rj/r. The exponent a of the exponent is plotted here as 
a function of St at R\ = 83, and for the values of e = 5 
(very low gravity) and eo = 0.4 (moderate gravity). The 
value of the exponent increases sharply as a function of 
St up to St « 1, where it starts to saturate and decrease 
slightly. The qualitative aspect of the dependence of a 
as function of St does not depend on the precise value 
of the Reynolds number in the range of R\ studied. We 
find that at values of St > 0.1, the value of the exponent 
is somewhat higher than the one found by Falkovich et 
al, 2002, see Fig. 4. This difference can be attributed 
to the fact that the exponent there was computed by 
studying contraction along the fluid trajectory, which in 
the limit St — > differs very little from the particles tra- 
jectories. Quantitative differences remain even for values 
of St as small as St R* 0.1. At moderate values of St, 
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FIG. 3: The blow up frequency as a function of the Stokes 
numbers for several Reynolds numbers: R\ — 45 (dot-dashed 
line), 83 (dashed line) and 105 (full line) and at eo = 5 and 
e = 0.4. 
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FIG. 4: The St-dependence of the exponents a, obtained by 
fitting (n 2 ) r by a power law dependence as a function of r : 
(n 2 ) r oc (rj/r) a , at a very small value of gravity (eo = 5, upper 
curve) and at a moderate value of gravity (eo = 0.4, lower 
curve). The value of the exponent obtained by Falkovich and 
Pumir (2004) at small values of St, is shown by the dashed 
line. 



the exponent is larger when gravity is small, as expected 
(Falkovich et al 2002). At larger values of the Stokes 
number (St > 2), the two exponents obtained with a 
very low gravity (e = 5) and with a moderate gravity 
(eo = 0.4) become very close to one another. 

Finally, Fig. 5 shows the continuous contribution to the 
collision rates, normalized by 327ra 3 r^-, as a function of 
the Stokes numbers. The continuous contribution starts 
at a low value at St — > 0, close to the value predicted by 
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FIG. 5: The Stdependence of the continuous component 
of the collision rate, normalized by the particle size, a and 
by the Kolmogorov time scale, tk, at small (eo = 5, upper 
curve) and moderate (eo = 0.4, lower curves) value of gravity. 
The horizontal dashed line corresponds to the Saffman and 
Turner formula. The upper graph corresponds to a Reynolds 
number of Rx = 45, the lower graph to a Reynolds number 
of R\ = 105. The collision rate peaks at a value St ss 1. 



the Saffman and Turner formula (indicated by the hori- 
zontal dashed line) , and increases to a maximum value at 
St w 1., before decreasing very slowly. Increasing gravity 
(decreasing eo) tends to decrease the collision rate. Over 
the range of parameters studied here, it was found that 
the continuous part of the collision rate increases when 
the Reynolds number increases. 

The sling contribution to the collision rate, shown in 
Fig. 6 at the value of the Reynolds number R\ = 105, 
starts from essentially zero at very small values of the 
Stokes number (the probability of having a sling effect 
is practically zero at St <C 1). Again, similarly to what 
has been observed for the continuous contribution to the 
collision term, it increases to a maximum at St ~ 0.8. 
The phenomenological description of the sling collision 
rate used in this work is not expected to hold at values 
of the Stokes numbers larger than ~ 1. In fact, our ap- 
proach is based on the implicit assumption that particles 




FIG. 6: The Sf-dependence of the sling component of the col- 
lision rate, normalized by the particle size, a, and by the Kol- 
mogorov time, tk, at small (eo = 5, upper curve) and mod- 
erate (eo = 0.4, lower curves) value of gravity, at a Reynolds 
number of R\ — 105. 



can be described by an essentially smooth hydrodynamic 
representation. This assumption becomes questionable 
as soon as St > 1. For this reason, only the part of 
the curve corresponding to values of St < 1.0 has been 
shown. Fig. 5 is meant to show the main trend, at moder- 
ate Stokes numbers (the formula used to define this term, 
EqOHis defined up to a constant). These data do not al- 
low to test the simple approximation St^ 1 ' 2 exp(— A) St) 
for the sling contribution suggested by Wilkinson et al 
(2006) for the case without gravity. 

Our method to estimate collision rates, although based 
on procedures which can be formally completely justified, 
is very indirect. It is therefore appropriate to compare 
our estimates for the collision rates with the results ob- 
tained by Franklin et al., 2005, by using direct numeri- 
cal simulations, and by estimating in a straightforward 
manner the collision rate among particles. We ran a sim- 
ulation at a value of R\ = 45 (run 3a) comparable to 
the run 3 of Franklin et al, at R\ = 48. The values of 
their 'Geometric collision rates', estimated for two par- 
ticles' sizes : a — lOfim (T\o) and a = 20/im (r2o), see 
their Table 5, are to be compared directly to the value of 
Kcont and of K s u ng , computed in the present work. The 
values obtained with our method, after proper rescaling, 
are shown in TableQ] At the lowest value of St, the value 
of r and K cont are very close. Part of the difference can 
be explained by the fact that our Reynolds is slightly less 
than the one obtained by Franklin et al. No sling contri- 
bution is expected in this case. At the highest value of 
the Stokes number, one finds that the continuous part of 
the collision rate, K cont , underestimates the value found 
in Franklin et al. On the other hand, a significant sling 
effect is expected at this value of the Stokes number, so 
the difference can be interpreted as resulting from the 
sling contribution. 

In conclusion, we have studied the collision rates in- 
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TABLE I: Comparison between the numerical estimates of Franklin et al., 2005, and the present estimates. The subscript 10 
(respectively 20) refer to particles of size lOiim (respectively 20/xm). 



Run # R x a St e r K cont K s i ing 

3 48 10pm 0.08 0.21 1.0 x 10~ 6 

3a 45 0.08 0.2 8.5 x 10~ 7 

3 48 20Atm 0.32 0.21 5.9 x 10 -6 

3a 45 0.2 2.2 x 10~ 6 1 x 10~ 6 



duced by turbulent air motion. The method used in this 
work is essentially lagrangian. We follow particles ad- 
vected in the flow, compute directly the flux of incom- 
ing particles (continuous contribution) and estimate the 
number of collisions occuring in the aftermath of a 'sling' 
effect. 

The ratio of the collision rate to the Saffman- Turner 
formula is found to increase significantly from 1 to ~ 10 
when the Stokes number increases from St ~ to St « 1 . 
The increase becomes more pronounced as the Reynolds 
number becomes larger. 

In the range of Reynolds number studied here, sling 
contributions are negligible at very small Stokes numbers: 
their probability goes as exp(— A/St) as a function of St, 
with a coefficient A of order 1. In practice, they become 



significant for Stokes numbers St > 0.20. 

The actual collision rates computed in this work are 
consistent, at R\ w 45, with the recent results obtained 
by Franklin et al. (2005). In particular, our results al- 
low us to disentangle the contributions due to the sling 
events, which we find to be quite significant for particles 
of size a — 20fim. 

This work should help clarify the origin of the enhance- 
ment of the collision rates of inertial particles due to tur- 
bulence, and also ultimately, to devise a parametrization 
of this collision rate. 
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